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A  procedure  for  adjusting  the  gains  in  an  a filter  used  in  tracking  air  targets  by  search 
radars  is  given  for  the  case  in  which  the  track  updates  appear  randomly  in  time.  The  filter 
gains  are  given  by  a  =  1  -  e'Z&oT  and  p  =  1  +  e-2*:;joT  -  2  e~$uoT  cos  wc’T  where  £,  ujq» 
and  ua  are  constants  and  T  U  the  randomly  vaiying  time  between  updates.  Using  this  gain 
adjustment  procedure,  we  found  that  the  tracking  errors  are  smaller  than  when  the  gpins  a 
and  p  are  held  constant  for  tracks  which  are  randomly  updated. 
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ABSTRACT 


A  procedure  for  adjusting  the  gum  in  an  a-0  £iter  used  in  track¬ 
ing  air  targets  by  search  radars  is  jpven  for  the  case  in  which  the 
track  updates  appear  randomly  in  time.  The  filter  gains  are  given  by 
a  =  1  -•  e~^uoT  and  0  =  1  ♦  e~2iuo'r  ~  2  e'&'i7  cos  uidT  where  £, 
wo,  and  are  constants  and  T  is  the  randomly  varying  time  be¬ 
tween  updates.  Using  this  gain  adjustment  procedure,  we  found  that 
the  tracking  errors  swe  smaller  than  when  the  gains  or  and  ?  are  *  dd 
constant  i  >r  tracks  which  are  randomly  updated. 


Manuscript  submitted  August  21, 1973. 


GAIN  ADJUSTMENT  OF  AN  ALPHA-BETA  FILTER 
WITH  RANDOM  UPDATES 


INTRODUCTION 

Search  radars  sometimes  track  air  targets  by  use  of  an  a-3  filter  (!}.  The  filter 
computes  the  target’s  velocity  from  the  measured  pcstion,  smoothes  both  the  p«  itior, 
and  velocity,  and  finally  predicts  the  position  the  target  will  have  at  the  nest  looL  of  the 
radar.  In  most  of  these  cases  a  uniform  update  tune  can  be  assumed  and  such  a  system 
can  be  analyzed  w:th  standard  techniques  involving  sampled-data  systems  (2-4}. 

Let  us  consider  merging  the  target  reports  from  a  number  of  radars.  The  radars  are 
assumed  to  rotate  at  different  speech,  and  the  detection  capability  of  any  given  radar  on 
a  target  depends  on  many  environmental  and  radar  factors.  For  these  reasons  the  track 
updates  for  a  given  target  might  be  thought  of  as  appearing  randomly  in  time.  Under 
some  circumstances  a  phase d-array  »adar  could  also  be  modeled  with  a  random  update 
time.  In  this  report  we  analyze  the  a-3  filter  by  using  the  postulated  random  updates, 
and  we  consider  a  means  of  adjusting  the  filter’s  gain  so  as  to  improve  its  system  response. 


REVIEW  OF  THE  a-0  FILTER 

The  review  of  the  a  —0  filter  with  constant  update  times  includes  an  examination  of 
the  transfer  function,  its  properties,  and  the  response  of  tlie  system  in  both  the  mean 
and  variance.  We  begin  by  defining  the  a-3  filter  (1): 


-  xp(k) 

+  o(xm(fe)  -  Jtp(fe)}  . 

(1) 

1»  +  ( 0fT)\xm{k )  -  xp(k)\  , 

(2) 

Xp(k  +  1)  -  Xflk) 

*  ut(k)T, 

(3) 

where 


x,(k)  =  smoothed  position, 
v,(k)  -  smoothed  velocity, 
xp(k)  -  predicted  position, 
x„,ik)  -  measured  position. 


T  =  sampling  period  (constant  until  specified  otherwise), 
a,  p  -  system  gains. 
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W.»  obtain  a  convenient  form  of  the  filter  equation?  by  substituting  (3)  into  (1)  and  (2), 
obtaining 


Xt(k  -  1) 

+ 

a 

•Mk  -  1) 

0IT 

[*.{*>1 


xp(*  +  l)  =  (IT)  ; 


Applying  the  z  transform  (2, 3)  to  (4)  and  (5),  we  find  that  the  transfer  functions  of  the 
system  are 


Hx  -  xt(z)jxm(z)  = 


Hu  =  vt(z)lxm(z)  = 


az  | z  +  (£  -  cr)/a] 
z 2  -  2(2-  a-  0)  +  (1  -  a) ' 

{g/T)c(r-4) _ 

z2  -  z(2-o~0)  -  (1-a)’ 

,  ..  ,  .  (a  +  0)ziz  -  (a/(or  +  0)]| 


By  setting  z  =  eA*,r)  we  find  .’.he  frequency  responses  for  a  typical  system  {Fig.  1).  The 
smoothed  position  is  obtained  by  passing  the  measured  position  through  a  low-pass  filter, 
and  the  smoothed  velocity  is  obtained  by  differentiating  the  measured  position.  The 
sampler  itself  acts  as  a  low-pass  filter,  and  any  excitation  whose  frequency  rang'  is  above 
1/7*  is  simply  folded  into  the  frequency  range  from  zero  to  1  IT.  The  normally  used 
values  of  a  and  0  are  shown  in  Fig.  2.  This  is  obtained  by  examining  the  pole  and  zero 
locations  of  the  transfer  functions  (6).  (7),  and  (8).  The  transfer  functions  (6),  {7),  (8) 
are  placed  into  standard  notation  for  a  second-order  system: 


H,.,  = 


z 2  -  2z  e~^w°T  cos  i *>dT  +  c~^°r 


Equating  terms  in  the  denominator  of  (3)  with  the  denominators  of  <6)s  (7).  and  (5),  we 
obtain 

o=l-  €'2^>°T  (10) 

and 

0  =  1  ♦  e'2^^  -  2e'^°r  cos  u)dT ,  (11) 

where  J  b  the  damping  coefficient  of  the  second-order  system,  Wq  is  the  natural  frequency, 
and  -.-:a  is  the  damped  natural  frequency.  Equations  (10)  and  {11}  will  be  used  later. 

Before  observing  the  system  response,  we  consider  incorporating  a  eonstar.-.-raie  high¬ 
speed  sampler  in  the  system  and  defining  as 


PARAMETER  0 
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Fit-  1— Frequency  response  of  an 
e-$  filter,  which  i»  given  by  £cjv. 
(6),  (7),  *nd  (S)  with  a  -  */«~ 
for  a  »  0.629  and  d  -  0.579 


F%,  2— PenniMthle  nio-*  of  a  and  f? 


BEK  a  CANTHELL 


*n(£)  =  u{tk)  +  ie(;*). 


with  w(tk)  assumed  to  be  »  whi .e  stationary  Gaussian  noise  having  zero  mean  and  having 
variance  gJ*  representing  measi jement  error  of  the  radar  and  with  uf  f *)  assumed  to  be 
sample  from  a  deterministic  target  trajectory.  A  block  diagram  of  the  filter  is  shown  in 
Fig.  3.  We  define  the  matrices  in  Eqs.  (4)  and  (5)  as 


AvT)  = 


(l-o)  (1-a.T 

-sir  a  -  p) 


.r-i. 


B<D=j  • 

m)  *  |  • 

jH 

cm  =  [i  idtj , 


•<**»  S*i 


xiu  _ _ 


JOTW-n 


U  *m  Li  rf!  iY  J  rCONSSWUT- 

i  II  ^  j  ImCM-SPES 


i  I  ■  j  \.KSGM-5PEE0Sawsi£f> ) 

F:s*.  3— An  0-0  fitter incorporating  a  cootinl-ritr  high-speed  sampler 

where  A  is  the  time  between  samples  of  the  highspeed  sampler.  We  then  write  the  filter 
equations  (4)  and  (5)  as 

Xfk)  =  A{T)X(k~  1)  +  B(D  [><(**)  +  U7(f*))  (13) 


xp(kT+i±)  =  C(t)X(k),  (14) 

where  i  =  1,  2,  3,  . . 77 A,  with  T/A  being  an  integer. 

The  response  of  the  filter  is  obtained  in  terms  of  means  and  covariarces.  The  equa¬ 
tions  describing  the  covariances  are  of  the  form  (4|. 


where 


P(fe  +  1)  =  A(77  P(Jt)  A’(77  ♦  B[T)  oj  Km  . 


~PIiXJk)  Fx^(k) 

l)  =  * 

Fp x  \kj  iff  (k) 

m-  *  *  *  *  _ 


in  which 


hut**.  ■iMifiHfiifu  ii  mu  . . mi  i  mi  iinttiiiiitnt(Miiiiliriii  n  iiiiiiftftiihinBiftiurBiitim**— J“,‘1 
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r«c.  i— Jiortnafixed  milnct  of  the  pr«u=t«d  poetkae  **. 
nrrjjJr  '.niiiu  i  it  a  facetiae  of  <?■  ied  p 


at  uniform  instsits  in  time  tccrstant  71,  the  response  of  tee  ;"-?ter  ■**{*),  *r><i 

xp<k7  ♦  T)  are  sinusoids  of  the  same  frequency  but  different  amplitudes  and  phases  at 
the  sampiLty  instants.  This  result  b  well  known  J2.3  j,  and  is  «^2v  obtained  tq?  re- 
cursreely  sorring  the  filter  equations  using  the  sampled  values  uUth  The  mtrssampb 
ripple  is  obtained  by  usng  the  constant-rate  h^h-speed  sampler  :  =  1.  2,  ....  7/a.  A 
system  response  including  the  aitrasampte  rip  pic  is  shown  in  Fig.  6.  Thb  filter  was  nsis- 
designed  purposely  so  that  the  error  could  be  shown  easily. 

The  various  ways  of  designing  the  filter  task  ally  attempt  to  reject  the  influence  of 
Es(f)  while  maintaining  Xp{k)  as  dose  as  possible  to  s*(t*}  1 1.6-10].  We  now  consider 
Use  filter’s  operation  when  it  b  randomly  updated. 
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GAIN  ADJUSTMENT  WITH  RANDOM  SAMPLING 

We  begin  by  ecmpuiina  the  means  and  covariances  of  the  filter  under  random  sam¬ 
pling,  considering  only  the  excitation  a>(tt)  for  a  fixed-gain  ssd  a  \=aub;*>gair,  a-£  filter. 
We  then  compute  the  me®  and  covariances  of  the  o--?  filter  considering  only  the  ©tdta- 
tiors  effi)  and  X?0>.  The  total  mean  and  covariance  response  can  be  found  by  super- 
unposing  tise  two  solutions,  A  procedure  for  showmg  that  superposition  e  pmaasfl^  is 
shown  in  Appendix  A. 


Fised-Csin  a-p  Filter 


The  a-p  fQter  with  random  updates  can  be  placed  in  the  form 

X{«1  =  AfJAjXCe-  I)  4  S(7ilfi.-{;it  +  w|ti>l 


f$5 


Xpir^.j)  -  ClTj.tlXliJ, 


where  Tk  »  a  random  variable  repeeatiling  the  toe  between  th*  Ik  -  Iitr  svi  Irth 
smnple,  -  Tj  ♦  T\  *  ...  *  TV  h  the  time  of  oectsTenee  oi  the  eth  sample,  »d  Air4h 
Bf  TV),  and  CiTi)  *re  matrices  whose  dements  depend  or.  7t.  The  joim  probability  dm- 
sHf  fi~)  of  the  process  i* 


flT-mr,)  ...  lX(0)i  /ImUjII 


fi* 


■«•'»>! /|afr,i,  aft,) 


Random  ssrspies  »t  (rim  a  conUnusus  determhsttbc  ail)  m*  ect  independent. 

The  time  ini«rral  T,  between  is  assumed  to  be  Kkmtxclir  Juthbuteu  for  afl  ; 

from  the  uniform  distribution  shown  at  fig.  The  subscript  on  Tk  wf 1  he  droned  for 
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E-  » 


t  d*« 

SAMPLING  PERIOD  ? 


Fig.  7—  Probability  density  of  a  sampling  interval  T 


The  covariances  of  the  system  states  an.  _*xt  considered  Considering  the  response 
only  due  to  w(tk),  we  find  X(fe)  =  0,  where  the  bar  denotes  expected  value.  As  shown 
in  Appendix  A,  the  recursive  equation 


P(k)  =  A(7)P(ft-l)A'(7)  +  8(7)  oJb’{T) 


(26) 


can  be  used  to  find  the  covariances  P(ft)  due  to  '-u(tk)-  The  covariance  equations  for  the 
a-/3  filter  are 


Px,x,(k) 

(1  "  o)2 

2(1  -  a)2T 

(1  -  a)2r2 

Px*(k-  1) 

Px,v,(k) 

- 

-0(1  -  a)/T  (1  -  or)(l  -  2)3) 

('l  -  flr)(l  -  0)  T 

Px,vt  (k~l) 

L 

(8  IT)2 

-20(1  -PUT 

(1-0)2 

Pv..A*- 1) 

L_  ~  -J 

a2 

+ 

a0/T 

(/?/T)2 

PXpxp(tk+ 1)  «  Px,x,ik)  +  2 TPXt.Jt(k)  +  T*PVtVt(k) 
Since  a  and  0  are  constant,  the  coefficients  in  Eqs.  (27)  and  (28)  are 


(27) 


(28) 


(Hd)  J  7  d7  =  (d/2)  +  e  , 

(29) 

T2  =  (d2/3)  +  de  +  e2 , 

(30) 

Iff  *  [Sn  (d  +  e)  -  2n  e]/d. 

(31) 

l/T2  =  l/(e(7  +  e)l  . 

(32) 

As  the  results  of  using  the  coefficients  (29)  through  (32)  in  (27.  and  (28),  the  normalized 
variance  of  the  predicted  position 
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°02  *'  ^xpxp(lk+l )l°Ji  (33) 

is  shown  in  Fig.  8  as  a  function  of  e  for  steady-state  conditions.  These  resuits  were  also 
computed  using  Monte  Carlo  techniques.  Observing  Fig.  3.  we  find  that  o02  increases 
rapidly  for  small  values  of  e.  This  can  be  explained  as  follows:  In  a  uniform  update 
system,  if  the  update  time  becomes  short,  tire  variance  in  the  velocity  Increases  rapidly. 
However,  this  effect  is  directly  canceled  in  the  variance  of  the  predicted  jiosition,  because 
one  needs  to  predict  only  over  this  same  small  interval  of  time.  In  the  random  update 
system  the  rime  between  updates  can  first  be  short,  creating  a  large  error  in  velocity. 

This  can  then  be  followed  by  a  long  time  interval  in  which  the  target’s  position  must  be 
estimated  using  the  poor  velocity  estimate.  Therefore  o02  increases  under  these  condi¬ 
tions.  A  method  of  avoiding  these  large  errors  is  next  considered. 


«  (SECONDS) 

Fig.  8— Normalijed  predicted  positicr.  variance  at* 
function  of  e  wi  'it  d  *=  6  seconds  /or  a  fixed-gain 
a-0  fvter 


Variable-Grin  a~0  Fflter 

The  method  of  avoiding  large  errors  in  the  random  update  system  adjusts  the  system 
gains  a  and  0  according  to  (10)  and  where  T  is  the  random  update  time.  The 
covariance  equations  considering  only  w{ as  ar.  excitation  is  again  given  by  Eqs.  (27) 
and  (28).  The  coefficients  are  computed  by  numerical  integration  with  the  use  of  the 
uniform  probability  density  given  in  Fig.  7.  The  normalized  variance  of  the  predicted 
position  (33)  is  shown  in  Fig.  9  as  a  function  of  e.  We  find  that  ot2  is  not  a  strong 
function  of  e,  as  was  the  case  in  Fig.  i.  This  can  be  explained  as  follows:  When  the 
time  interval  between  samples  becomes  short,  cr  and  0  approach  0,  thus  smoothing  the 
data  heavily  and  negating  the  rise  in  velocity  errors.  As  the  time  between  samples  be¬ 
comes  long,  a  and  0  approach  1  and  in  effect  no  smoothing  is  used.  The  gain  adjustments 
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Fig.  B— Nortiwlired  predicted  politic  •  variance  as  a 
function  of  €  with  {  «  0.4  and  d  *  6  seconds  for  a 
variable-gain  a-j}  filter 


(10)  and  (11)  appear  to  maintain  a  rather  constant  o02  under  the  random  sampling  dis¬ 
tribution  proposed. 

We  now  investigate  the  effects  of  the  excitation  u(fj-)  and  X(0)  on  the  system  re¬ 
sponse  using  the  gains  (10)  and  (11).  We  will  find  the  effect  of  u(tfc)  and  X(0)  on  the 
randomly  updated  filter  by  using  simulation  procedures.  (Because  E tlu(fy)  -■  u(f/)] 

(u(f,-)  -  u ((,•)](  is  not  G  for  all  i  *  j ,  an  analytic  formulation  appears  to  be  very  difficult) 
An  approach  to  the  problem  would  be  to  find  the  mean  and  covariance  of  the  predicted 
position  at  the  hth  sample.  However  the  feth  sample  appears  randomly  in  time.  A  more 
meaningful  calculation  would  be  to  find  these  quantities  at  a  given  instant  of  time.  Con¬ 
sider  the  constant-rate  high-speed  sampler  defined  in  Fig.  3.  Undo'  any  randomly  sampled 
excitation  the  filter  response  can  be  found  at  a  given  sample  of  the  high-speed  constant- 
rate  sampler  which  corresponds  to  a  given  instant  of  time.  The  quantity  A  is  made  smaller 
than  e  such  tliat  the  effect  of  each  random  sample  can  be  easily  seen.  The  adjustable-gain 
filter  is  excited  by  taking  random  samples,  distributed  ««  -I.own  in  Fig.  7,  from  a  sinusoid, 
given  by  Eq.  (22).  Over  many  trials  the  mean  and  variance  of  xp  is  computed  at  each 
switch  closing  of  the  constant-rate  high-speed  sampler.  The  results  are  shown  in  Figs.  10 
and  11.  Even  though  u(<)  is  a  deterministic  process,  under  random  sampling  the  filter 
output  is  a  random  variable.  One  finds  that  if  the  filter  has  a  sufficient  bandwidth,  the 
niter  with  adjustable  gains  (10)  and  (11)  follows  u(?)  quite  well  (Fig.  11). 

The  reason  xp  follows  u(t)  fairly  well  with  the  adjustable-gain  filter  may  be  explained 
as  follows:  As  the  time  between  samples  becomes  short,  the  system  smoothes  heavily, 
which  counteracts  the  increase  in  bandwidth  due  to  the  rapid  sampling.  Conversely,  as 
the  time  between  samples  becomes  large,  the  system  lightly  smoothes,  counteracting  the 
sluggishness  induced  by  the  long  time  between  samples.  In  fact  the  system  attempts  to 
maintain  a  constant  bandwidth,  and  if  this  bandwidth  is  wider  than  the  frequency  content 
of  the  signal,  then  xp  follows  u(t)  fairly  closeiy.  This  system  behaves  like  a  continuous 
system  of  fixed  bandwidth  being  excited  with  inputs  appearing  randomly  in  time. 
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Fig.  10— Mean  (and  mean  p!us  and  minu*  •  .  jndird  deviation)  of  xp 
vs  time  for  ^  •  0.4,  w0  “  0.1  raa/s,  d  *  8  *,  and  e  "  0.1 1 


Fig,  1 1  —  Mean  (and  mean  plus  and  minus  a  standard  deviation)  ol 
xp  n  time  for  $  *  0.4.  u>0  =  0.314  rad  It,  o'  «  6  s.  and  {  »  0.1  s 
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V> 

The  total  filter  response  in  both  the  mean  arid  covariance  can  be  obtained  by  super¬ 
imposing  the  effects  of  and  (X(0).  u(f)l-  The  filter  design  would  adjust  £  and  uig 
such  that  xp  followed  u(t)  as  dosely  as  {vossible  while  minimizing  Pxpxp>  which  is  com¬ 
posed  of  the  variance  Pxp-p  due  to  tv{t)  plus  the  variance  Pxpxp  due  “10. 


SUMMARY 

The  constant- coefficient  a-ti  filter  when  randomly  updated  was  found  io  have  large 
errors  in  the  variance  under  certain  conditions,  namely,  a  short  time  between  updates 
followed  by  a  long  time  between  ^updates.  To  circumvent  this  problem,  variable  gains 
a  -  1  _  e_2$wo T  and  £  =  l  +  t.*2^oT  _  2  e~h^oT  u>dT,  where  £,  <Oo,  and  are 
co  ns  tar'  and  T  is  the  randomly  varying  time  between  updates,  were  postulated.  Using 
these  gains,  one  found  that  the  variance  in  the  predicted  position  remained  at  reasonably 
low  values  under  the  same  conditions.  Also  it  was  shown  that  the  variable- gain  filter’s 
response  to  a  given  trajectory  could  achieve  a  reasonable  small  error. 

Although  the  gain  adjustments  were  simply  postulated  and  shown  to  work  well  with 
the  examoles  cited,  no  exact  justification  for  their  use  was  given.  However  arguments 
were  given  as  to  why  the  system  seemed  to  work  well. 
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Apoendix  A 

FORMULATIONS  OF  MEAN  ,  ND  COVARIANCE  EQUATIONS 

The  purpose  of  this  appendix  is  to  cutliiw  a  procedure  for  showing  that  superposi¬ 
tion  can  be  used  in  computing  the  means  and  covariances  of  (23)  and  to  outline  a  method 
of  obtaining  (26).  We  begin  by  writing  the  system  equations  (23)  using  a  convenient 
form. 

X(fe)  =  AX(*-1)+  Blu(ffc)  +  io(tfc)J  .  (Al> 

where 

A  =  random  variable  which  is  identically  distributed  and  independent  from 
sample  to  sample. 

B  -  random  variable  which  is  identically  distributed  and  independent  frorr. 
sample  to  sample. 

u(hfe)  =  random  samples  from  an  arbitrary  function, 

X(0)  =  initial  condition. 

iqtir)  -  zero-mean  white  Gaussian  noise  with  variance  ou,2. 

Recursively  solving  (Al),  we  can  obtain  the  solution  in  the  form 

X(h)  =  g|  A,  B,  X(0) .  w(tfc) ,  u>(f*)]  .  (A 2) 

The  mean  M(ft)  ant'  covariance  P(«)  of  X(k)  are  computed  using  the  probability  density 
(25).  The  means  M’(&)  and  M’(ft)  and  the  covariances  P^fr)  and  P “(fc)  are  computed  in 
the  same  manner  by  first  using  X(0)  and  u(ffc)  as  excitations  and  then  using  iv(.‘*)  as  an 
excitation.  Superpostion  is  then  shown  to  hold  by  noting 

M{k)  =  M'(fe)  *  M”(fe),  (A3) 

P(k)  -  P'(Ar)  t  P"{k)  .  ;A4) 

The  computation  is  straightforward  but  quite  lengthy. 

We  compute  (26)  by  noting  the  mean  value  of  X(fc>  excited  with  only  is  0 

and  forming  the  covariance: 

X(fejKlfej  «  A  X(k  -  li  X'(k  -  1?  A7  +  1 fuTtkU^iW 

+  AX(D-  l)tt(y  f  +  B  it»(fjj)  X(fe  -  1 )  A'  .  (A5» 
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where  the  bar  denotes  expected  value  and  ifie  prime  denotes  transpose  Beo.ise  of  the 
independence. assumed,  (AS)  become 


X(*)  X'(ft)  =  A  X(*t  -  1)  X'(fc  -  1)  A*  +  B  w(tk )  w(tk)  & 


+  A  X(fe  -  1)  w(tk)  3*  +  B  &{tk)  X(h  -  1)  A'  .  (A6) 


Defining  P(fc)  =  X(fc)  X*(k)  and  noting  X(fe  -  I)  ic(t*)  =  0,  we  obtain 


P(*)  =  A  P(A  -  1)  A'  +  B  o2  B‘ .  (A7) 


